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Abstract. Monte-Carlo simulations for a resonance scan of the charmonium-like state X(3872) at PANDA are performed. 
Final state radiation hadronic background reactions are taken into account. The signal reconstruction uses a realistic pattern 
recognition (track finder and track fitter) and electron/pion discrimination. 
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THE X(3872) STATE 

The X(3872) is a charmonium(-like) state which was first observed in its decay X(3872)— hF/y/TT+TC - in B decays [1] [2] 
and inclusive production in pp collisions [3] [4] . Meanwhile more decays such as final states D°D°* or J/xj/co have been 
found. In particular, the observation of its radiative decay into Jlyy allows the assignment of positive charge parity 
C=+l. The interesting observation of isospin violation in the decay X(3872)^//y/p(^7T + 7C~) raised the question, if it 
might be a charmonium state at all. As its mass is within Am< 1 MeV of the D°D * threshold, it might be interpreted as 
an S- wave molecular state [6] [7]. Many properties of the X(3872) are being studied at ongoing experiments. However, 
the width is unknown and probably can only be determined at PANDA. The current upper limit is T<2.3 MeV [1], 
whereas PANDA might be able to set an upper limit in the order of a few hundred keV. The quantitative investigation 
of this upper limit on the width by detailed Monte-Carlo (MC) simulations incl. background is topic of this paper. 

CHARMONIUM PRODUCTION AT PANDA 

The PANDA (Anti-Proton Annihilation at Darmstadt) experiment will be located at the future FAIR (Facility for 
Antiproton and Ion Research) facility [8] at GSI in Darmstadt, Germany. The primary synchrotron will have a 
circumference of 1.084 km, a proton beam with 2>Q<Ei, eam <90 GeV and an intensity <2xl0 13 /s. Antiprotons are 
produced by a secondary target and then stored and cooled in the HESR (High Energy Storage Ring). There are 
two HESR modes. In the high intensity mode (high resolution mode), there are 10 11 (10 10 ) stored antiprotons, and 
stochastic cooling (electron cooling) will provide a beam momentum resolution of Ap/p>l-10~ 4 (Ap/p>4-10~ 5 ). 
PANDA will investigate pp collisions with an antiproton beam momentum p/, eam <15 GeV/c, corresponding to 
,/s<5.5 GeV. In e + e~ collisions, only states with quantum numbers J PC =1 can be produced directly due to the 
intermediate virtual photon. As the X(3872) has positive C parity, thus, it can not be produced directly in e + e~. 
However, in pp collisions two or three gluons are produced in the annihilation, and states with any quantum number 
can be formed. 

At the PANDA interaction region, three different types of targets are foreseen, namely (a) hydrogen pellets with a 
thickness of ~25 jim and a falling velocity of >60 m/s, (b) a cluster jet target and (c) fixed nuclear targets such as 
Be, C, Si or Al. For the pellet target, a peak luminosity of Jz?=2x 10 32 cm~ 2 s _1 is expected. If adjusted to a single 
resonance, such as e.g. the J/y/ this would correspond to about 2x 10 9 produced J/y/ per year. 

A resonance scan with the cooled beam could provide a measurement of the X(3872) width. This technique was 
pioneered by the Fermilab experiments E760 and E835. With a beam momentum resolution of Ap/p=2-10~ 4 , 



corresponding to a ^ resolution of 0.5 MeV (FWHM), the measurement of the width of the J/y with 
r(7/i//)=99±12(stat.)±6(syst.) keV and the width of the \j/ with r(y/)=306±36(stat.)±16(syst.) keV could be 
performed [9]. 

The requirements for the PANDA experiment are quite high: capability for high rate detection and data acquisition for 
<2x 10 7 interactions/s, c±4n solid angle coverage (incl. a high resolution forward detector, as PANDA is a fixed target 
experiment), secondary vertex detection of e.g. D mesons with a resolution a<100 jUm, charged track momentum 
resolution Ap/p~l%, charged particle identification for e ± , fj, ± , n ± , K ± , p ± using e.g. Cherenkov detectors with 
internal reflection, and electromagnetic calorimetry in the range 10 MeV<£y<5 GeV. Details about the experiment can 
be found elsewhere [10]. For the down below results, MC simulations were performed with all PANDA subdetectors 
implemented, but the reconstruction used only the following PANDA subdetectors: 

• a Time Projection Chamber (TPC) with 135 padrows (corresponding to < 135 hits for a charged particle track) 
and 135,169 pads of 2x2 mm 2 area, 

• a Micro Vertex Detector (MVD) with 120 silicon pixel modules (lOOx 100 jUm 2 pixel size and in total 10 7 readout 
channels), and 400 silicon strip modules (~0.5 m 2 active area and 7x 10 4 readout channels), and 

• an Electromagnetic Calorimeter (EMC) with ~ 17,200 crystals of PbWC>4, a radiation hard scintillation material 
with a fast decay constant of ~6 ns. The thickness corresponds to ~28Xo radiation lengths. 

The MC simulation was performed with the PandaRoot [11] simulation, digitization, reconstruction and analysis 
framework. It consists of ~43,000 geometry volumes (incl. details such as e.g. the beampipe for the pellet target) 
and >400,000 lines of C++ code. Transport engines are Geant3, Geant4, and Fluka. Event Generators are EvtGen, 
DPM, PYTHIA, and UrQMD. References can be found elsewhere [12]. PandaRoot is used on >15 Linux platforms. 
Recent improvements, in particular since the PANDA Physics Report [13] are in particular: (a) for the first time 
X(3872) simulations, not investigated before as a physics channel for PANDA, (b) usage of detailed field maps of a 
homogenous B z =2 T in the central region and dipole field of 2 Tm in the forward region (with flux effects in the iron 
of the PANDA muon subdetector and interference of the fields taken into accout), (c) usage of a realistic track finder 
and track fitter based upon a conformal map technique, and (d) simulation of final state radiation [14]. In about ~30% 
of all //y/^e + e~ decays an additional photon is radiated. 

The particle indentification uses the ratio Elp as a variable for discrimination of charged pions and leptons. E is the 
deposited shower energy in the EMC, p is the reconstructed track momentum by the track finder and track fitter using 
the MVD and the TPC. For leptons Elp is ~1. For charged pions, which generate hadronic showers, Elp can take any 
value 0<E/p<l. By applying an Elp cut of 0.8 <Elp< 1.2 for electrons, the combinatorial background is reduced by 
>90%. 

ESTIMATED RATES FOR X(3872) FORMATION AT PANDA 

For resonant X(3872) formation, the required beam momentum is pt, eam =6.99100 GeV/c. Our baseline assumption 
is a peak cross section of <7x(3872)=50 nb, in the same order of magnitude as e.g. the cross section of y/ production 

[9]. Note, that, if the X(3872) is a molecular D°D°* state, there are estimates [15] that the cross section may increase 

up to <443 nb. The ratio of branching fractions of the X(3872) decays into the final states D°Z5 *:7/v/^ + ^ _ :7/i^7 
is assumed to be 9:1:0. All other decays are assumed to have a zero branching fraction. The branching fraction for 
J/\jf— >e + e~ and ji + pL~ is ~6% each. The reconstruction efficiency of ~50% is dominated by the track reconstruction 
efficiency for the low momentum charged pions. All these factors lead to a reconstructable cross section of 250 pb. 
For resonance scans, the HESR will be operated in the high resolution mode with Ap/p=10~ 5 , which corresponds to a 
luminosity of «5f=2x 10 31 cm~ 2 s~'. Assuming an accelerator duty factor of 50%, this leads to an integrated luminosity 
of ^=0.86 pb-Vday. For a resonance scan with 20 energy points and 2 days/point, this would correspond to a yield 
of ~215 events of pp^X(3E72)^J/y/n + n^ per day at peak. 

BACKGROUND 

The main background is meson production in processes such as pp^n + n~n + n~ with two misidentified charged 
pions (as leptons) in the EMC. The cross section for this process is 50 jub [16], compared to an estimated signal cross 
section of 50 nb. For comparison, the total pp cross section is ~70 mb. It is important to investigate the shape of the 



background e.g. in the 2-particle mass spectrum. For this purpose, a DPM (dual parton model) event generator [17] 
[18] [19] was used. In fact, a varying background shape in the region of the J/\j/ was found and fitted with a first order 
Chebyshev polynome. The fit was performed differently for each beam momentum in the resonance scan. 



DETERMINATION OF A RESONANCE WIDTH BY RESONANCE SCAN 

The Breit-Wigner cross section for the formation and subsequent decay of a cc resonance R of spin J, mass Mr and 
total width Fr formed in the reaction pp^R is 

(27+1) Anjhcf T 2 R BR(pp -> *) x BR(R -» f) 

(2S+l)(2S+l)(E 2 n -4(m p c 2 ) 2 ) X (E cm — Mrc 2 ) 2 + F r /4- " W 

where S is the spin of the (anti-)proton. 

a(E cm ) = [ a BW (E')G(E' - E cm )dE' 
Jo 

is a convolution of a Breit-Wigner term for the resonance and the function G for the beam resolution. 
The area under the resonance peak is given by 
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which importantly is independent of the form of G(E). C pea \ is the cross section at E cm =M R c 2 given by 

(27+1) 16nh 2 BR(pp —>•/?) x BR(R — > /) 
a " eak ~ (25+l)(25+l) (M R -4m 2 )c 2 . ( } 

By measuring A using a fit to the excitation function and inserting <J pea k into Eq. 3, the resonance width Fr can be 
determined. 



PRELIMINARY RESULTS 

MC simulations for pp— >X(3872)— >J/\j/n + 7t^ with background were performed for 9 different beam momenta in the 
resonance region. Fig. 1 shows the e + e~ invariant mass with PID applied and background subtracted. The final state 
is exclusive, i.e. no other particle is produced except for photons from final state radiation (visible as radiative tail at 
masses smaller than the Jly mass). Each data point corresponds to 2 days of data taking. A Gaussian fit was used 
to extract the signal yield. This Jly yield is equal to the tagged X(3872) yield. Then data points were fitted with an 
excitation function. The preliminary results are as follows: 

• The fitted width of the excitation function is ~20% larger than the width which was used as input in the MC 
simulation, due to the analysis technique such as the background subtraction and the reconstruction technique 
(e.g. the track fitter). Note, that the beam resolution is not the limit for the width measurement. Even with a beam 
momentum resolution of Apb ea m— 0.5 MeV, the E760 experiment was able to determine the factor of ~5 smaller 
width of the Jly [9]. The difference between generated and reconstructed width would only be a systematic error 
on the width measurement. 

• The width of the excitation function is still <100 keV, if the assumed resolution of the anti-proton beam is 
A/?//?=210 -5 (HESR high resolution mode). However, in this analysis all distributions were assumed to have a 
Gaussian shape, while the correct resonance shape would be given by a Breit-Wigner parametrisation and the 
shape of the function G(E) in Eq. 2 might be non-Gaussian. In addition, for our analysis, we approximated the 

function V^ = \j ^ m p + ^ m p \J Pbeam +™%by a linear function in the small range of the resonance region. As long 

as the beam resolution is accurately known, our results indicate that the width of the X(3872) could be determined 
at PANDA, if it is larger than r^>100 keV and if the DPM generator describes the background shape correctly. 

These results are still preliminary, as the unfolding (integral equation, Eq. 2) is ongoing work. 
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FIGURE 1. MC simulation of a resonance scan of the X(3872) at PANDA, with 9 scan points of different anti-proton beam 
momenta and data taking of 2 days for each point. The plots show the background subtracted, tagged J/y/ signal from the 
X(3872)^y/y/7T + 7T~ decay. Only statistical errors are shown. 
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